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Reversible modifications of linear dispersion - graphene between boron nitride 

monolayers 
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Electronic properties of the graphene layer sandwiched between two hexagonal boron nitride sheets 
have been studied using the first principles calculations and the minimal tight-binding model. It is 
shown that for the ABC-stacked structure in the absence of external field the bands are linear in the 
vicinity of the Dirac points as in the case of single layer graphene. For certain atomic configuration 
the electric field effect allows opening of a bandgap of over 230 meV. We believe that this mechanism 
of energy gap tuning could significantly improve the characteristics of graphene-based field effect 
transistors and pave the way for novel electronic applications. 
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I. INTRODUCTION 

Graphene, the thinnest electronic material with amaz- 
ing properties^ is a promising candidate for next- 
generation electronics. Its cone-shaped valence and con- 
duction bands touch each other at the corners of the Bril- 
louin zone. The large area graphene is then semimetallic 
and not suitable for use in digital logic devices due to 
the fact that such channels cannot be switched off. How- 
ever, the band structure of graphene can be modified and 
many approaches have been proposed to break the sym- 
metry between sublattices responsible for massless char- 
acter of Dirac fermions and the lack of the bandgap. 2 
Theoretical investigations have demonstrated that a nar- 
row bandgap can be created by placing graphene sheet 
on the lattice-matched materials, such as hexagonal 
boron nitride (h-BN). 3 An alternative route to induce 
the bandgap is functionalization of graphene with suit- 
able elements, as fluorine (fluorination^) and hydrogen 
(hydrogenetion^), which transforms graphene into an in- 
sulator with bandgap of about 3.5 eV. This process 
is reversible, but it takes hours to regenerate pristine 
graphene from these materials. 

On the other hand, energy gap tuning by electric field 
effect in bilayer Bernal stacked graphene has been widely 
studied. e,T Recently, a bilayer graphene transistor with a 
sizable on/off current ratio has been fabricated 8 showing 
a great potential of graphene-based electronics. One of 
the main difficulties of producing bilayer graphene tran- 
sistors is disorder. For example, graphene devices on 
Si02 substrate are highly disordered due to the scat- 
tering from charged surface states and impurities which 
limits the carrier mobility and suppresses the graphene 
unique properties. The measured moblities 10000-15000 
cm 2 V -1 s -1 are much smaller than those reported for 
suspended graphene. However, the very recent experi- 
mental studies have revealed^ that graphene transferred 
onto h-BN substrate is free of wrinkles, distortions and 
surface charge traps. This leads to the construction of 
graphene devices with very high mobility close to 60000 
cm 2 V" 1 s- 1 . 



The conventional semiconductors exhibit a general 
trend that with increasing bandgap, the electron mobility 
decreases. 9 This behavior has been also predicted for car- 
bon nanotubes and graphene nanoribbons. In both above 
mentioned cases, the valence and conduction bands be- 
come progressively parabolic (rather than linear) when 
the larger bandgap opens. It causes the increase in the 
effective masses of charge carriers and, in consequence, 
the decrease in their mobility. In the unbiased bilayer 
graphene the low energy bands are parabolic, not cone- 
shaped, which lead to a finite value of the effective mass. 
If the linear dispersion of pristine graphene was pre- 
served in the tunable system, it would offer new possibil- 
ities for electronics based on the properties of massless 
Dirac fermions. In particular, the high carrier mobility 
in graphene at room temperature would profit without 
the loss of the possibility to switch devices off. 

Our recent theoretical studies^ on graphene/h-BN 
systems have revealed that the energy gap tuning in 
the range starting from to 130 meV can be achieved 
in graphene/h-BN bilayer system. However, the linear 
dispersion is not exactly preserved - the mass of Dirac 
fermions is small, but cannot be reduced to zero. 

In this paper we report on the possibility of creating 
and tuning the bandgap in trilayer systems consisting of 
graphene monolayer sandwiched between two h-BN hon- 
eycombs (see Fig. [TJ. It seems that in ABC stacked tri- 
layer h-BN/graphene/h-BN system the presence of two 
identical boron-carbon or nitrogen-carbon dimers in the 
unit cell does not induce the high asymmetry between 
sublattices, which guarantees the conical dispersion ex- 
actly preserved in the K points. When the electric field 
is applied perpendicularly to the layers, the bands are 
shifted and deformed which leads to the bandgap open- 
ing. 



II. GEOMETRY OF TRILAYERS 

We have studied a stability of graphene/h-BN AB- 
stacked bilayers and h-BN/graphene/h-BN trilayers with 




(b) 




FIG. 1. (Color online) (a) Lattice structure of h-BN/graphene/h-BN ABC-stacked trilayer system, darker and lighter balls 
denote atoms of the same type in h-BN layers (b) orientations of graphene and h-BN layers with respect to each other. 



ABA and ABC stacking orders by means of density func- 
tional theory (DFT) calculations^. In case of the bilayer 
there are two inequivalent orientations of the graphene 
with respect to the h-BN lattice. In the lowest-energy 
configuration carbon atoms Bi and A2 are placed directly 
above the boron atom {A\ ) and on top of the empty cen- 
ter of the hexagon in the bottom layer, respectively. In 
the second configuration, the boron and nitrogen atoms 
are exchanged, thus the unit cell contains a nitrogen- 
carbon (N-C) dimer. 

The relative position of h-BN and graphene layers al- 
lows for two different orientations of the third layer in 
the sandwich. The same orientation of both h-BN lay- 
ers gives ABA Bernal stacking configuration typical for 
graphite. In the ABC (rhombohedral) stacking (see Fig. 
[T] (b)) each pair of nearest-neighboring layers forms an 
AB-stacked bilayer with the upper B sublattice directly 
on top of the lower A sublattice and the upper A above 
the center of the hexagon below J^ The optimized dis- 
tances between the layers are the same for ABA and 
ABC stacking orders. The difference in total energies 
is less than 0.1 meV which is within a computational er- 
ror. Moreover, it seems that the particular configuration 
would be determined by the substrate, as in the case of 
the graphene trilayers on SiCJ^ In case of two inequiva- 
lent configurations containing N-C or B-C dimers, there 
is a larger difference in total energies and interlayer sep- 
arations, but both of them are stable. The energy differ- 
ence is then over 30 meV preferring the latter configura- 
tion, while the distances between the layers are equal to 
3.44 and 3.22 A, respectively. 



III. TIGHT BINDING APPROACH 

We have analyzed the electronic properties of h- 
BN/graphene/h-BN trilayers using the nearest-neighbor 



tight-binding model. We include only the interlayer hop- 
ping between atoms which are on top of each other in 
adjacent layers, thus the interaction between the sheets 
takes place through half of the atoms in each plane. This 
method has been previously applied to study trilayer 
graphene and allowed for simple and qualitatively cor- 
rect analysis of this system, with and without the applied 
external fields. 15 

First, we remind the hamiltonian describing the bot- 
tom bilayer in the system. We include the parameter U 
of the interlayer asymmetry induced by an external field, 
which is uniform and perpendicular to the layers. In case 
of the AB-stacked bilayer we construct the hamiltonian 
in the basis with atomic components $Ai , $Bi , &A 2 > ^b 2 
(only 2p z orbitals of respective atoms are taken into ac- 
count, that is &Ai is B 2p z orbital with Bloch factor, $Bj 
N 2p z , $a 2 j &b 2 are 2p z orbitals of carbon atoms with 
Bloch factors)ii: 
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where we fixed the zero energy as 2p level of carbon 
atoms in the crystal. 

As in the Slonczewski-Weiss-McClure model of the 
AB bulk graphite, Jq (70) and 71 describe the nearest- 
neighbor intralayer and interlayer hopping parameters 
respectively, e a and e* 3 denote onsite energies of in- 
equivalent atoms in the h-BN honeycomb and f(k) — 

^2i=i elkSi > where 8 stands for the in-plane nearest- 
neighbor vectors. 

We generalize this minimal model on the case of the 
ABC trilayer including &a 3 and $b 3 components and 
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consider the case when the lower layer is at potential 
U/2, the upper layer is at potential -U/2 and the middle 
is at zero potential: 



/ 

where we neglected all hopping parameters between non- 
adjacent layers^. 

It is worth to note that the above hamiltonian re- 
sembles that of ABC graphite^ and graphene trilayer—. 
However, the system contains only one graphene layer 
and its properties are more similar to graphene monolay- 
ers rather than to trilayers. 

In the corners of the Brillouin zone, |/| =0, thus the 
eigenvalues can be easily determined 
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One can observe that in the absence of the external 
field (U — 0) there are three double roots of the secular 
equation. It means that there is no energy gap between 
valence and conduction bands, which is clearly visible 
in Fig. [5] (a). DFT calculations have shown that the 
degeneracy is in part broken by next-nearest neighbor 
interlayer hopping between sites A\ and B3 that lie on 
the same vertical line^S, but it does not affect the highest 
valence and the lowest conduction bands. 

Moreover, it is easy to show by the methods of the 
perturbation theory that there are two pairs of bands, 
which are linear in the vicinity of the Dirac points: 



Ei, a {k) = £1,2 ± 



7o|/(fc)l(7Jh) 



(6) 



since |/(fc)| is linear there. It has to be stressed that the 
system is half-filled and the Fermi level coincides with 
the touching point of the cones. 

The valence and conduction bands are linear in the 
wide range of energies (±1 eV with respect to the touch- 
ing point) with the value of Fermi velocity nearly equal 
to that of graphene. The external electric field (U 7^ 0) 
lifts the degeneracy in the K point leading to shifts and 




FIG. 2. TB band structure of h-BN/graphene/h-BN ABC- 
stacked trilayer with N-C dimers in the absence (a) and in 
the presence (b) of the external electric field perpendicular to 
the layers. 
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FIG. 3. (Color online) Variation of the energy gap as a func- 
tion of U in the range, where the gap is direct. The case of 
ABC-stacked trilayer containing N-C dimers is chosen. Red 
(dashed) line denotes the dependence for biased graphene bi- 
layer. For simplicity the same value of interlayer hopping 71 
is taken for both systems. 



the repulsion of bands (see Fig. [21(b)). The variation of 
the induced energy gap between valence and conduction 
bands is shown in Fig. [3] and compared to the bandgap 
opening in biased bilayer graphene. 



IV. NUMERICAL RESULTS AND DISCUSSION 

The ground state configuration containing B-C dimer 
does not offer the possibility of energy gap tuning in the 
sufficiently large range, which is related to the onsite en- 
ergies values. Among h-BN/graphene/h-BN sandwiched 
trilayers, the ABC-stacked system containing nitrogen- 
carbon dimers seems to have the largest potential for ap- 



plications. We have found the approximate values of the 
tight-binding parameters describing this system^ The 
onsite energies values are e a = —1.50 eV and e^ = 3.16 
eV representing nitrogen and boron 2p levels in the crys- 
tal lattice. Intralayer hopping parameters are equal to 
7o=2.28 eV and 7o=2.59 eV for h-BN and graphene 
layers, respectively. The value of an out-of-plane hop- 
ping has been estimated to 7i=0.25 eV. According to 
the predictions of TB model with this set of parame- 
ters, the largest energy gap that can be created between 
the valence and conduction bands is nearly 230 meV and 
strongly depends on the 71 value. DFT results show that 
the bandgap achieves 232.6 meV suggesting that the pa- 
rameters slightly change in biased system. Figure |4] il- 
lustrates DFT predictions for the band structure of the 
system with the effective value of U exceeding 3 eV, where 
a large bandgap is opened. 



V. PERSPECTIVES 

In summary, we have explored that in the ABC- 
stacked h-BN/graphene/h-BN trilayer one can observe 
the conical energy dispersion relation near the Dirac 
point, which is typical for pristine graphene. The electric 
field effect allows to create an energy gap up to about 230 
meV depending on the exact value of interlayer hopping 
parameter 71. 

We proposed a unique mechanism of bandgap tuning in 
graphene-based system, which, under certain conditions, 
could be an alternative for graphene bilayer on Si02- Fig- 
ure 3 illustrates that in the considered trilayer the gap 
opens slowly as compared to bilayer graphene. DFT re- 
sults show that the screening effects are of the same order 
as in the previously studied graphene/h-BN bilayer— as 
well as in the bilayer graphene^i, thus, in practice, the 
opening of sizable bandgap can be difficult. On the other 
hand, the disorder which is partly responsible for diffi- 
culties of producing bilayers transistors, in graphene/h- 
BN systems is significantly reduced indicating the very 
high quality of this substrate^ In the present work, we 
have studied the issue of the band structure modifica- 



tions. It is very important, but not the only factor that 
determines the characteristics of graphene FETs. We 
believe that the advantage of linear dispersion preserved 
in the unbiased trilayer can improve current on/off ra- 
tios or the speed of graphene devices^, however, a fur- 
ther theoretical analysis of the real systems and exper- 
imental verification have to be performed. Preliminary 
results on fabrication of graphene/h-BN combined struc- 
tures have been recently reported ) 10 ' 22 thus we hope that 
our work will encourage experimentalists to produce h- 
BN/graphene/h-BN stacks. 




FIG. 4. DFT band structure along TKM line in momentum 
space of trilayer h-BN/graphene/h-BN in the presence of ex- 
ternal electric field applied in the direction perpendicular to 
the layers. The indirect energy gap is 231.5 meV. 
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